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Preface : 

This paper is a transcript of a talk entitled "Black Holes - XI: 
'X-Ray Observations" given by one of us (E. B. ) at the "International 
Conference on X-Rays in Space" at the University of Calgary (Canada) 
August 14-21, 1974. Xt served as the companion lecture to one by 
Dr, R. Ruffini (Princeton) entitled "Black Holes - I: Theory". 

Although the results from several key experiments are reviewed this 
presentation emphasizes recent work by the Goddard X-ray astronomy 
group, extending the analysis of our data on Cygnus X-l beyond the 
rapid temporal variations reported earlier (Rothschild et al., 197*0. 
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WHAT IS SPECIAL ABGJT CYGNUS X-l? 

(BLACK HOLES '01 THEORY AMD OBSERVATION II: X-RAY OBSERVATIONS) 

Elihu Boldt, Stephen Holt 

Richard Rothschild* and Peter Serlemitsos 

Goddard Space Flight Center 
Greenbelt, Maryland 

ABSTRACT 

Of the eight X-ray sources now known which may be associated with 
binary stellar systems, Cygnus X-X is the most likely candidate for 
being a black hole. The X-ray evidence from several experiments is 
reviewed, with special emphasis on those characteristics which appear 
to distinguish Cygnus X-l from other compact X-ray emitting objects. 

Data are examined within the context of a model in which millisecond 
bursts (Rothschild et al,, 1974) are superposed upon shot-noise 
fluctuations (Terrell, 1972) arising from "events" of durations on the 
order of a second. Possible spectral-temporal correlations are inves- 
tigated, indicating new measurements that need to be made in future 
experiments. 

INTRODUCTION 

A more accurate title for this talk might be "What is Special About 
Cyg X-l?", For comparison in e> .his, Table 1 lists the eight 
X-ray sources currently discussed' as being members of binary stellar 
systems. The ordering is according to binary period, ranging from the 
4.8 hours for Cyg X-3 to times of about a week or more. All but two 

*NAS-NRC Resident Research Associate 
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of these exhibit pronounced occulintions. Of the remaining two Cyg X-3 
exhibits a large but smooth variation in intensity for the few keV band, 
with a well established periodicity (Parsignault er al, , 1972; Canizares 
et al, , 1973)* Since the spectrum of Cyg X-3 shows severe absorption 
effects (cf., Bleach et al, , 1972 ) at a few keV and the hard X-ray flux 
may not be periodically modulated (Baity et al., 19735 Ferettl et al., 
197 *+ ) , observed variations are likely to be associated with obscuration 
of the source, However, for Cygnus X-l, the other non-occulting 
source, there are no strong absorption effects apparent in the few keV 
region of the spectrum although there is some evidence for intensity 
variations with a 5.6 day period (Sanford ?.t al., 197*0. Hence, Cyg X-l 
may be unique in the sense that it affords us a relatively unobscured 
view of the accretion process close to the compact object where observ- 
able differences between a white dwarf, neutron star and black hole 
could be most pronounced (Ruffini, 197*+). 

Of course, our strong interest in Cyg X-l is motivated by its 
apparently high mass. If the binary interpretation is correct (cf., 
Bahcall et al. , 197*+) and if the distance to the optical member is not 
vastly overestimated (cf . , Margon et al. , 1973) , this mass is the only 
one listed in Table 1 which exceeds 3 solar masses, the upper limit 
for a neutron star (Leach and Ruffini, 1973). Therefore, the interpre- 
tation of Cyg X-l as a collapsed object implies a black hole. Effects 
which could be supporting evidence for this interpretation are: 

1) X-ray emission driven by an accretion process 

2) Absence of sustained oscillations in the emission (cf . , Leach 
and Ruffini, 1973) 
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3 ) Chaotic fluctuations (cf . , Shakura and Sunyaev, 1973) down 
to the millisecond region for a Schwarz child metric black 
hole and to the sub-millisecond region for an extreme Kerr 
black hole. 

The 2-20 keV spectra for the sources listed in Table 1 are char- 
acterized as either "flat" or "steep" except for the case of Cyg X-l. 

This exception is related to the situation that the X-ray emission from 
Cyg X-l sometimes exhibits a highly dynamic spectrum in the 2-20 keV 
band. This is probably indicative of an accretion process with special 
characteristics of the kind to be explored by further studying sp< sctral 
correlations with radio emission and optical emission such as the \ 4686 
line from He II (cf,, Hutchings ot al, , 1973) in the accreting gas. 

Except for the pulsars Her X-l and Cen X-3, sustained oscillations 
(i,e., pulsar activity) seem to be absent for all sources listed in 
Table 1, to the limits indicated. Although an observation of Cir X-l 
in 19 69 by Margon et al. (l97l) detected a 1.5 Hz oscillation present 
during their 8 second exposure, subsequent observations with UHURU 
(Jones et al, , 1974) and an MIT rocket-borne experiment (Spada et al., 
1974) did not see any evidence for a pulsar. However, sporadic flaring 
on time scales down to about a second were evident, such that low duty 
cycle quasi-periodicities can not be ruled out. For the casr of Cyg X-l, 
rapid sporadic behavior is most pronounced, manifested by millisecond 
bursts superposed on large "shot* noise" variations that are rich in 
quasi-periodic behavior at frequencies in the vicinity of 1 Hz. These 
temporal effects for Cyg X-l will be exhibited in detail later. 
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Finally, in comparing the sources listed in Table 1, we should not 
omit the fact that Cyg X-l is different from the others in one obvious 
way that could give us an observational bias in that its peak intensity 
is the highest of all, However, I should also note that our observation 
of millisecond bursts from Cyg X-l (Rothschild et al. , 1974) occurred 
during a rocket-borne experiment when the average intensities observed 
for Cyg X-3 and Cyg X-l were within about 20$ of each other, 

II, SPECTRUM 

When Cyg X-l was first observed with UHURU in late 1970 it showed 
a much steeper spectrum than had been seen with rocket-borne experiments 
during a span of several years preceeding 1970. As shown in Figure 1, 
it was not only unusually steep but a number of sightings on the same 
day (December 21, 1970) gave clear evidence that this spectrum was highly 
variable (Schreier et al,, 1971). Except for a Goddard observation in 
September 1970 (Bleach et al., 1972) and an NRL observation in June 1964 
(Bowyer et al., 19^5 ),« all measurements prior to UHURU were essentially 
consistent with a non-variable flat spectrum in the few keV region. 

This relatively stable situation seems to have been restored in March 
1971 with the sudden onset of detectable radio emission associated with 
Cyg X-l. The remarkable correlation between the X-ray and radio trans- 
itions is exhibited in Figure 2, which shows how the hard X-ray emis- 
sion (> 10 keV) increased and lower energy X-ray emission (2-6 keV) 
decreased when the radio source became evident (Tananbaum et al., 1972). 
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The average spectral situation for Cyg X-i observed vdth UHURU 
prior to the radio onset way he compared with that obtained almost a 
year later by examination of the two curves exhibited in Figure 3. 

All measurements since 1971 that I am aware of indicate a spectrum 
comparable to the flat one (energy spectral index of about l/2) in 
this figure. However, a soft X-ray measurement by Stevens et al. 

(1972) in October 1971 may be interpreted (Garmire, 197*0 in terms 
of a steep component decreased in intensity from that Eteep component 
shown in Figure 3, so that it becomes recognizable only below' about 
a keV. In this way, the amount of interstellar X-ray absorption 
needed to match the soft part of the spectrum becomes high enough to be 
compatible with that recently inferred by Ryter et al, (1974) on the 
basis of interstellar reddening. Since Cyg X-l may have undergone a 
radio intensity variation close to the time of the soft X-ray observation 
in October 1971 (Hjellming, 1973) 5 possible correlations of any steep 
X-ray spectral component with relatively small radio transitions sub- 
sequent to the large one in March 1971 need to be pursued further. 

The spectral history of Cyg X-l is summarized in Table 2 in terms 
of measurements that bracket most closely in time the known pronounced 
spectral transitions. The high interisity of Cyg X-l in the 2-10 keV 
band observed by NRL in 1964 (Bowyer et al., 1965) was several times 
higher than that observed in subsequent rocket-borne experiments and 
comparable to that seen with UHURU just before the transition in March 
1971. Some time between the Livermore measurement in November 1968 
(Mac Gregor et al., 1970) and the Goddard measurement (Bleach et al. , 
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1972) in September 1970, the few keV spectrum began to steepen. For 
energies above 20 keV (i.e. , the part of the spectrum accessible from 
balloon-borne observations) two of the sharpest transitions were a 
three fold decrease from May to June 1967 (Overbeck and Tanaribaum, 1963) 
and a six fold decrease from June to July 1969 (Matteson, 1971). It 
was not until the UHURU observations during March-April 1971 that a 
hard X-ray increase was observed, on a comparable time scale, and where 
we learned of an apparent anti-correlation of the spectral transitions 
above and below about 10 keV. 

Ill . 5-6 DAY VARIATIONS ? 

Although there seem to be some variations in the X-ray emission 
from Cyg X-l associated with the binary period, the situation is apparently 
more complicated than it is for the other well established binaries. 

Figure 4 summarizes observations of Cyg X-l made with the University 
College (London) experiment on 0 AO- Copernicus (Mason et al. , 1974), 
giving the exposure according to the radial velocity phase of the 
spectroscopic binary HDE 226868 (zero phase corresponds to the secondary 
on the far side of the primary) as well as the time measured in Julian 
days. For eight passes through or near zero phase, there were five 
occasions when the intensity showed a dip lasting an hour or two. For 
the first four of these events, the measured spectrum indicated increased 
absorption; the fifth dip associated with zero-phase may be indicative 
of the chance of an associated coincidence with random intensity fluc- 
tuations since this "event" showed no absorption. At 0S0-7 experiment 
(Li and Clark, 1974) examined Cyg X-l at zero-phase on four occasions * 
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and observed an absorption event once (also shown on Figure 4); the 
characteristic energy for unit optical depth changed from about 0.6 keV 
to about 1.9 keV during the dip. So, the score seems to be five clear 
absorption events out of twelve that would be expected for a strictly 
regular phenomenon. Furthermore, the 0A0 observations suggest that these 
events are not synchronized to zero-phase and may be progressively 
shifting in phase. In summary, these events are recurrent but not 
regular, and it appears that while the binary motion in itself is a 
necessary condition it is rot sufficient for their occurrence. Mason 
et al. (1974) interpret this behavior da terms of a model whereby the 
observed X-ray emission is intercepted by a variable stream of matter 

that bridges the primary star and the X-ray source. 

( 

The last continuous Copernicus observation of Cyg X-l shorn in 
Figure 4 corresponds to an interval of 6.5 days in November 1973 which 
extended through two zero-phase transitions. The overall intensity 
pattern during this observation suggests a broad maximum centered at 
phase 0 . 5 , with a maximum to minimum intensity in the ratio 4/3. How- 
ever, considering a series of zero-phase observations spread through- 
out the preceeding year, the maximum to minimum intensity is also in 
the ratio 4/3. I conclude, then, that the suggestion of a 5.6 day 
periodic intensity variation for Cyg X-l is in need of more supporting 
evidence. 

VI. FAST FLUCTUATIONS 

The difference between Cyg X-l and other X-ray emitting members 
of binary systems as regards short-term variations is illustrated in 
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Figure 5. The figure exhibits samples (15 seconds each) of the count- 
rate profile for Cyg X-l, Gyg X-3 and Her X-l obtained during the 
same rocket-borne Goddard experiment on October 4, 1973 (Rothschild 
et al. , 19?4) , Even without folding, the raw Her X-l data clearly 
show the 1,24 a period pulsar, To within the statistical uncertainties, 

Gyg X-3 shows no signal variations on any of the time scales examined. 

In sharp contrast, Gyg X-l exhibits rapid, large and apparently chaotic 
fluctuations . 

The first high temporal resolution study of Gyg X-l after the spectral 
transition of March 1971 was made on May 1, 1971 with a rocket-borne 
experiment by Rappaport et al. (1971). This observation lasted some- 
what over a minute and involved 1 millisecond resolution. The count 
rate profile ,“y r this exposure is shown in Figure 6 and is qualitatively 
similar to that obtained over two years later in the Goddard experiment. 

The low frequency portion of the power density spectrum for these MIT 
data is shown in Figure 7, Although the portion of the spectrum above 
~ 1 Hz gave no indication of significant deviations from what would be 
expected from white noise and rocket motion could moke some contribu- 
tion below ~ 0,2 Hz the remaining region (i.e. between 0.2 Hz and 1 Hz) 
exhibits so much power that this frequency band alone corresponds to 
about a 3C$ modulation of the total signal. 

In terms of apparent periodicities, the detailed temporal behavior 
of Cyg X-l seems to change from one observation to the next. The 
discovery of rapid pulsations for Cyg X-l by Oda et al. (1971) suggested 
there might be important frequency components up to ~ l4 Hz during 
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various 20 second duration UHURU exposures. An 8 second rocket obser- 
vation by Holt et al, (l97l) indicated frequency components up to ~ 5 Hz. 
Subsequent longer duration observations with UHURU (Sehreier et al, , 1971) 
indicate no persistent periodicity but the short term presence of many 
frequencies ranging from ~ 0,1 Hz to ~ 3 Hz and persisting for times on 
the order of a few seconds. Figure 8 is a sonogram (Oda et al,, 1972) 
of the MXT rocket data (Figure 5) and shows how various frequency Compon- 
ents up to ~ 5 Hz (200ms) come and go during that 70 second exposure. 

Is this quasi-periodic behavior genuine, indicative of some 
oscillatory physical process, or is it merely a manifestation of some 
other type of underlying variation? 

V. SHOT-NOISE MODEL 

Terrell (1972) examined the UHUBU data of Oda et al. (1971 ) in 
terms of the fluctuations arising from the superposition of randomly 
occurring pulses, (i.e. shot-noise). As shown in Figure 9? the 
computer simulations of the count rate profile based upon a random 
sample of shot-noise at an average rate of 400 sec" 1 for pulses of 
0.25 second duration each easily replicates the qualitative appearance 
of the UHURU data, 

A key aspect of Terrell* s analysis may be seen in Figure 10, where 
he exhibits th'e power density spectrum for the computer simulated count 
rate profile of Figure 9j significant power in excess of that expected 
for counting statistics alone is clearly evident. The important result 
is that shot noise adds greatly to power density fluctuations expected 
from just counting statistics and can actually masquerade as quasi- 
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periodic behavior in the frequency domain of interest. An elaboration 
of this effect may be seen in Figure 11 and Figure 12, Figure 11 shows 
the power density spectra for three IMMi passes at Cyg X-l and illustrates 
how, although the dominant frequencies fluctuate from pass to pass, the 
region of maximum powetf remains bounded at frequencies less than a few 
Hz, Figure 12 shows the power density spectra for five of Terrell’s 
computer generated random simulations of data using the same shot noise 
parameters (i.e,, 400 pulses sec. -1 , 0,25 sec pulse duration) and the 
temporal resolution and normalization of the UHCMJ data. Here again, 
the peak frequencies of significant power fluctuate from trial to trial, 
but are bounded within a few Hz. The solid curve in Figure 12 gives 
the average power density spectrum expected for the composite of shot 
noise and counting statistics and shows how the power increases at the 
low frequency end. The enhanced fluctuations in power corresponding 
to this expected average increase appear as significant power when erron- 
eously compared to the power expected from counting statistics alone 
(i.e., given by the dotted line in Figure 12). This then is an effective 
mechanism for generating the apparently quasi-periodic behavior seen 
for Cyg X-l. 

An integral measure of shot noise may be obtained by considering 
the variance in the observed count, as summarized in Figure 13. The 
shot noise is here parameterized by X (average pulse rate) and t ( tem- 
poral duration of the basic rectangular pulse). The data are represented 
by noting the counts (w) in each bin (width of AT) and constructing 
the variance C(N - N) 3 ] and the average (n). As indicated in Figure 13 s 
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if there is no shot noise the variance per unit bin width is expected 
to be just the average rate (i.e,, N/AT), independent of AT. However, 
if the data are dominated by shot noise with t«AT, the variance per 
unit bin width exceeds the average rate by a factor which depends upon 
\ but is again independent of AT. On the other hand, for shot noise 
dominated data with t»AT the variance per unit bin width increases 
from that for the case of no shot noise by a factor which goes linearly 
with AT. 

We have calculated the variance per unit bin width for the data 
of our 50 second exposure to Cyg X-l as well as the data for Cyg X-3 
obtained in the same flight (October 4, 1973) and the results are plotted 
in Figure 14. Except for statistical fluctuations that vary from the 
few percent expected for AT < 0.1 sec to the tens of percent expected 
for AT > 1 sec., the variance per unit bin width does indeed match the 
overall average rate for Cyg X-3. On the other hand, for Cyg X-l the 
deviation from the mean rate increases significantly with AT (for 
AT < 1 sec.) in a manner consistent with the initial linear behavior 
expected for shot noise. Also, we note that as AT > 1 sec. the variance 
per unit bin width begins to flatten, with variations dominated by 
statistical fluctuations. The behavior seen here may be characterized 
by \ = 19 sec. -1 average rate for shot noise pulses of duration T = l/2 
sec. This corresponds to an average of about 10 shot noise pulses over- 
lapping at any given time, a small enough number to give us an immediate 
qualitative picture for the large amplitude fluctuations in the Cyg X-l 
count rate we observed on time scales comparable to l/2 sec. 
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In order to analyse our Cyg X-l data for quasi-per iodic charac- 
teristics, the 5C second exposure was divided into 10 second intervals 
and the principal frequencies were identified, The power density spectra 
for the hand 10 Hz - 1.6k Hz gave no indications of significant deviations 
from what would be expected from white noise. However, as shown in 
Figure 15, below 10 Hz there was much significant power for Cyg X-l, 
especially in the vicinity of 1 Hz. For the threshold power considered 
here for identifying a principal frequency white noise gives a probability 
of 0.05 for obtaining an entry anywhere on this plot whereas the number 
of entries actually exhibited in Figure 15 is over a thousand fold that 
expectation. A similar analysis for Cyg X-3 exhibits no entries at this 
power level and gave no indication of significant deviations from what 
would be expected from white noise anywhere in the power density spectrum. 

With the 0.1 Hz resolution corresponding to 10 second intervals 
of data, several of the principal frequencies shown in Figure 15 for 
Cyg X-l appear to persist for consecutive intervals. However, up to 
frequencies of about 2 Hz there is so much power that about half of the 
possible frequencies qualify as principal frequencies and the random 
chance of repitition becomes appreciable. The histogram at the right 
hand margin of Figure 15 shows the rich concentration of principal 
frequencies in the band up to about 2 Hz and again illustrates the form- 
idable shot- noise associated background against which one must contrast 
any frequency in this range that may be a candidate for some genuinely 
periodic process on time scales ^ 10 seconds. 
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Adding our results for data obtained in 1973 to those tabulated 
by Terrell (1972) for observations Made through 1971 j the shot noise 
parameters for Cyg X-l can b<t traced through its recent history, This 
is summarized in Table 3> which gives the pulse rate (X), pulse width 
(t), and spectral state measured, for Cyg X-l on six separate occasions 
during the period 1967-1973, The spectral state is characterized here 
as being "high" or "low" and refers to the emission below 10 keV , as 
discussed before. The most remarkable aspect of the evolution so traced 
is the relative constancy of the characteristic pulse width, at a few 
tenths of a second, for pulse rates that differ by as much as a factor 
of forty. On the other hand, the pulse rate appears to be strongly 
related to the spectral state, where "high" corresponds to higher pulse 
rates. 

VI. PULSED OSCILLATIONS 

Since .the power spectrum for shot-noise dominated data is equal 
to the power spectrum for the underlying characteristic shot-noise pulse, 
we should be able to determine if these pulses carry higher frequency 
signals than represented by the much lower frequency components at a 
few Hz that characterize the pulse envelope. In other words, we should 
be able to tell if we are in fact detecting randomly pulsed oscillations. 
Since these shot-noise pulses last only a fraction of a second, the 
inherent resolution in locating carrier frequencies is a few Hz, and 
"tuning in" involves locating principal frequencies on a time interval 
comparable to the pulse width. 
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An analysis of the hind just’ described was performed on the Cyg X-l 
data of Goddard Flight 13.010 and the results are summarized in Figure 
16. Principal frequencies above 10 Hz are here indicated for consecutive 
intervals of 0.3277 s duration, comparable to the underlying shot-noise 
pulse width, for all of the Cyg X-l exposure. The frequency grid has 
a resolution of 3 Hz. For this analysis, the threshold power for defining 
a principal frequency is such that we expect white noise to exceed this 
threshold anywhere in the spectrum with a probability of 0.5 for each 
of the temporal intervals. First, we note that the total number of entries 
on this plot is consistent with that expected for white noise and that 
the population of such principal frequencies is statistically consistent 
with a flat distribution. Even so, the occasional repetition of a prin- 
cipal frequency (+ 3 Hz) in adjacent intervals could still indicate that 
a significant "carrier" frequency is sometimes present. However, this 
sort of repetition in consecutive intervals occurred just once and can 
not be taken as a statistically significant indication of such an effect. 

Figure 16 also shows the average count rate for each of the inter- 
vals used for this analysis and it may be of some interest to note that 
the peak intensity interval, containing three bursts of millisecond dur- 
ation (Rothschild et al., 1974), exhibited no principal frequencies 
above 10 Hz. Considering the entire exposure, the 99 $ confidence upper 
limit (cf. , Holt et al., 1973) to the fraction of the basic shot-noise 
pulse that may be modulated by frequencies in the band 10 Hz to 1.6 kHz 
is about 20$. 
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VII. TEMPORAL- SPECTRAL- IMMENSITY CORRELATIONS 

To examine possible short-term spectral-temporal effects for Cyg 
X-lj Brinkman eb al. {1974 ) examined UH11RU data in terms of an auto- 
correlation function. They separately considered two energy bands, 

2.1 - 5.1 keV and 5.1 - 12 keV, and the result for the lower energy 
band observed during a single 50 s pass is shown in Figure 17. The 
autocorrelation is evaluated at 0,192 s intervals of the time delay 
parameter. Three features to note are: 

1) a rapid decrease in the autocorrelation with the time delay 
parameter during the first three intervals (i.e. ^ 0.6 s) 

2) a slower less pronounced autocorrelation lasting about 
fifteen intervals (i.e. ~ 3 s) and 

3) possible quasi-periodic behavior for longer times. 

Comparing the high and low energy autocorrelation plots obtained 

from the composite of several sightings Brinkman et al. (197*0 noted 
an asymmetry exhibited in Figure 18. The short-term autocorrelation 
appears to persist about one time delay interval (i.e. ~ 0.2 s) longer 
for the low energy band than for the high energy band. 

We have examined data from Goddard Flight 13.010 in the same way, 
using energy bands and time delay intervals as close to those of 
Brinkman et al. (197*0 as readily feasible. The results of this analysis 
are shown in Figure 19 and are quite comparable to those of Brinkman 
et al, (1974), including the apparent difference noted between the low 
and high energy data. For comparison, we looked at the autocorrelation 
evaluated for the Cyg X-3 data in exactly the same way and obtained 
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the plots shown in Figure 20, a restilt to be expected for a flux dominated 
by statistical fluctuations. 

In order to gain some insight into the connection between the auto- 
correlation behavior and the shot-noise model, we computer simulated a 
count-rate profile based upon \ = 20 shot-noise pulses per second of 0.5 s 
duration each and used these simulated data for the evaluation of the 
autocorrelation, exactly as done previously for the real data. As 
shown in Figure 21, this model does a good job of simulating the short- 
term dominant autocorrelation present for the real data and shows us 
that the apparent quasi-periodicity in the autocorrelation first noted 
by Brinkman et al. (1974) is in fact just another artifact of the under- 
lying shot-noise. The simulation was repeated for shot-noise pulses of 
0.1 s duration and although the qualitative behavior is similar, a quan- 
titative comparison (i.e., correlation time and apparent frequency of 
quasi-periodicity) with that for the actual data (i.e,, Figure 19) shows 
that such a small pulse duration is clearly unacceptable. In neither 
simulation did we get any indication of the longer-term (~ 3 s) fall-off 
to the weaker portion of the autocorrelation observed in the Cyg X-l data. 

After comparing our results with those of Brinkman et al, (1974) 
using the same resolution for the increments of time delay, we repeated 
the evaluation of autocorrelation using an order of magnitude smaller 
resolution (i.e., 20.5 ms), and the results are shown in Figure 22. 

The fast and slow fall-off to the Cyg X-l autocorrelation are now 
clearly discernible as two components. The rapid one extrapolates to 
zero in 0.5 s whereas the slow one goes to zero at approximately 3 s. 
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As also shown In Figure 22 , we used this finer temporal resolution 
to again computer simulate the autocorrelation behavior. For an 
assumed shot-noise pulse width of 0,5 s we note that the autocorrelation 
goes to aero at 0.5 s and then oscillates with an average period of 1 s. 
Now that we see this simple situation we realize it should have been 
expected since* after all, the Fourier transform of the autocorrelation 
is the power density spectrum and, althougn the frequency spread for 
a 0.5 s duration pulse is 2 Hz, the dominant frequency component for 
such a pulse is in ^act Just 1 Hz. 

In order to see how critically tuned the high resolution auto- 
correlation is to the underlying shot-noise pulse width, we computer 
simulated two more cases, one for a pulse width of 0.1 s and the other 
for a pulse width of Q .98 s, and we can see (Figure 22 ) they are both 
clearly unacceptable, 

We have already noted that the autocorrelation function appears 
to be a function of energy; both the data of Brinkman et al. (197*+) 
and also our data from Flight 13.010 indicate that the principal portion 
of the autocorrelation persists somewhat longer (~ 0,2 s) for the low 
energy band (2,1 - 5*1 keV) then it does for the higher energy band 
(5.1 - 12 keV). To pursue this further we have also considered 
our data above 12 keV, even though the reduced count rate presents some 
quantum limitation difficulties in computation for the finest resolution 
(20.5 ms). An intermediate resolution (0.082 s) was used for comparing 
the autocorrelation at the highest energies (> 12 keV) observed for Cyg 
X-l with that for the two lower energy bands (see Figure 23) previously 
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considered, as well as with that evaluated for Cyg X-3, serving as a con- 
trol. The major result of this analysis is that the principal short- 
berm autocorrelation seems to have disappeared at the highest energies 
or, if there is a major short-term autocorrelation still present, it is 
shorter than the temporal resolution used here (i.e,, decays away in 
less than 0,082 s), 

In summary, essentially all the temporal variations for Gyg X-l 
on time scales from a fraction of a second to minutes and possible longer 
seem adequately described by shot-noise pulses of about 0,5 s duration, 
with the admixture of somewhat longer weaker components lasting up to 
a few seconds. Since the underlying shot- noise pulse width is relatively 
independent of source intensity we might expect that tracing the longer 
time scale variations of Cyg X-l In terms of \ could yield new evolution- 
ary information. That is, the long-time scale behavior should be charac- 
terized by the ratio of the variance per unit time to the square root of 
the mean rate as well as by the average rate Itself, Finally, since 
the shape of the basic shot-noise pulse itself appears to be a function 
of energy, the need for a more detailed spectral-temporal analysis is 
clearly indicated. As more data become available we should be able to 
reconstruct the spectral evolution of these individual events we have 
referred to as shot-noise pulses and thereby begin to understand their 
origin(s) and use them as probes of the otherwise inaccessible regions 
of closest approach to the compact- object at the core of Cyg X-l, 
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VIII . MILLISECOND BURSTS 

‘Hie characteristic pulses underlying the shot-noise behavior of 
the Cyg X-l fluctuations could themselves be the carriers of pronounced 
fluctuations on a smaller time scale (e.g,, millisecond bursts). This 
kind of sub- structure is suggested by two pieces of evidence. First, 
we recall that the M.I.T, rocket-borne experiment of 1971 May 1 (see 
Figure 6 ) had millisecond resolution; an analysis of that data by 
Oda et al. (1974) indicates that millisecond fluctuations are "bunched" 
to a degree that exceeds statistical expectations. Second, for the 
Goddard experiment (Rothschild eb al., 19?4 ) that involved a large 
enough detection area to identify 8 individually ree* gnizable millisecond 
bursts during a 50 s exposure 3 of them occurred within an interval of 
enhanced emission lasting about a third of a second (i.e., a time com- 
parable to a basic shot-noise pulse). 

As discussed earlier (see Section VI ), a power density spectrum 
for the Goddard data in the interval of enhanced emission containing 3 
of the 8 bursts identified does not yield any principal frequency (10 Us 
- 1,6 K Ha) and that implies an upper limit of about 20$ to the fraction 
of this enhanced emission which may be modulated by a periodic signal. 

To make a meaningful comparison with the observed "bur stlnes s", though, 
we need to obtain an estimate of the fraction of the emission in this 
pulse residing in millisecond bursts. In effect, this involves deter- 
mining how many bursts were actually missed using the rather stringent 
criterion for statistical significance that was applied to the entire 
exposure in identifying the 8 bursts already indicated. 
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To determine the "burst" fraction, we considered an interval of 
20 telemetry frames (i,e, , 0,4 s) encompassing the enhanced emission 

region where we found evidence for bunching and divided this into 320 
consecutive bins of 1,20 ms duration each, Noting the counts per bin 
as an independent variable, we obtained a histogram (Figure 24) of the 
observed population (cf , , Boldt et til,, 1971), Comparing this distri- 
bution with what would be expected from Poisson statistics (also shown 
in Figure 24), it is quite evident that the 4 bins with counts of 10 
or more each are not part of the average population, In fact, the 3 
bins of highest count are the ones previously identified as bursts (see 
Rothschild et al. , 1974) and all reside in one telemetry frame (i.e., 
within an interval of 20,48 ins,). If we eliminate this one telemetry 
frame from the 20 considered here and compare the resulting "purged" 
histogram with the now expected Poisson distribution (Figure 25 ) we 
see a much improved match. Furthermore, by examining Figure 25 in detail 
it becomes clear that we are left with a burst of 10 counts per bin and 
that there may still be anywhere from 0 to 8 more bursts with between 
6 to 8 counts per bin remaining, buried in the statistics. For the 4 
large unambiguous bursts alone, the burst fraction (i.e., burst counts/ 
total counts) for the 0.4 s unexpurgated intenral considered is 5$. If 
we add the additional 8 bursts which we might have missed we obtain an 
upper limit of 13 $ for this burst fraction. 

Now we are ready for the question of periodicity in connection 
with these bursts, The count rate profile encompassing the telemetry 
frame containing 3 bursts as well as adjacent frames is exhibited in 
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Figure 2 6 . We note that the time interval between the firct 2 bursts 
is 5 - 7 milliseconds, comparable to the period expected for the inner- 
most stable orbit around a 12 - l6 MS) Schwarzachild metric black-hole. 

For our present purpose, though, the importance is that we may now 
restrict our search for periodicities to frequencies within about a 
oO Hz band- width at about 170 Hz. And therefore, we may lower our 
previous estimate of the upper limit to the fraction of enhanced emis- 
sion modulated by a periodic carrier. However, the 3cj upper limit of 
l4$ so obtained is really no better than the 13$ upper limit to milli- 
second bursting (periodic or not) we obtained by simply extracting 
"would-be" bursts from the Poisson distribution of fluctuations. 

In summary, the fraction of a second "events" dominating the tem- 
poral fluctuations (i.e,, shot-noise) of Cyg X-l appear to be correlated 
with the bunching of millisecond bursts. However, the burst fraction 
for the one such event clearly indentified is not much greater than 
about 5$ and at most 13 - l4%. And we still can not say much about 
the question of periodicity except that wt seem to be on the threshold 
of experiments sensitive enough to answer this, perhaps even on the 
submillisecond time scale to be examined in our current rocket-borne 
experiment for Cyg X-l, 
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Figure 1 


Figure 2 


Figure 3 


Figure 4 


FIGURE CAPTIONS 

(a) Spectral density A versus spectral index for a large 
number of UHURU sightings of Cyg X-L on 1970 December 21. 
Dashed line, an approximate fit to the points. 

(b) Average spectrum of Cyg X-l on December 21, Solid 
lines, best-fit spectral indices to all the points and to 
the 10-20 keV data separately. Broken lines, typical bal- 
loon observation power-law fits. (This figure from Schreier 
et al., 1971.) 

The X-ray intensity of Cyg X-l in three energy ranges between 
2 and 20 keV and the radio fluxes at two frequencies are 
plotted against the day of 1970 (from Tananbaum et al., 1972), 

The average spectra for Cyg X-l from before and after the 
1971 March and April transition, as observed with UHURU 
(Tananbaum et al,, 1972). The spectra are plotted in pho- 
tons (cmP sec keV)" 1 versus energy, and various power law 
fits are labeled with the energy index. 

The observing periods of the seven 0A0 Copernicus observa- 
tions (Mason et al., 1974) are shown on a plot of radial 
velocity phase of the spectroscopic binary HDE 226868 versus 
time. The five Copernicus events are marked with the symbol 
X (the first four are absorption events). The OSO-7 M.I.T. 
observation of an absorption event (Li and Clark, 1974) is 
shown as an erro. bar corresponding to the estimated pre- 
cision in establishing the occurrence at zero phase. 
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Figure 5 


Figure 6 


Figure 7 


Figure 8 


Figure 9 


The £ount rate profiles* binned every 20,48 ms, for Her X-l, 
Cyg X-3 and Cyg X-l during 15-second exposure samples of 
each source obtained during the same rocket-borne experiment 
on 1973 October 4 (Rothschild et al., 1974). The inset 
shows the overall flight profile in counts per 2 second 
interval versus time after launch (in seconds). The arrows 
indicate where the 15-second data samples occurred relative 
to the full exposure. 

Cyg X-l counting-rate profile obtained during an M.I.T. 
rocket-borne experiment on 1971 May 1 (Rappaport et al. } 
1971). 

The low-frequency portion of the Cyg X-l power density 
spectrum for the data described in Figure 6 (Rappaport 
et al. , 1971) . 

A "Hissagram" (bottom) of the Cyg X-l data described in 
Figure 6 (reproduced at top) showing wave amplitudes as a 
function of wave period and time (i.e. a contoured sonogram) 
for 2,1 sec samples, obtained by Oda et al. (1972). 

A typical set of Cyg X-l counting-rate data from UHURU 
(bottom), due to Oda et al., (1971), together with computer- 
generated shot-noise data (top), due to Terrell (1972). 

The dashed line indicates a least-squares fit to the tri- 
angular trend of the observational data. 
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Figure 10 Raw, unsmoothed, power spectrum calculated from the shot- 
noise run shown in Figure 9 (top) , showing the large fluc- 
tuations at low frequencies typical of shot-noise spectra. 

The dashed line indicates the spectrum expected for counting 
statistics alone. (Terrell, 1972). 

Figure 11 Power spectra calculated from three sets of data due to 

Oda et al. (1971), normalized to 1.0 for counting statistics. 
(Terrell, 1972), 


Figure 12 Power spectra calculated from five sets of computer-generated 
shot-noise including counting fluctuations, normalized to 
1.0 for the part of the spectrum due to counting, indicated 
by the dashed line. The theoretical expectation is shown 
as a heavy line. (Terrell, 1972). 

Figure 13 A scheme for the integral measure of shot-noise parameters. 


Figure 14 Variance per unit bin width versus bin width (seconds) 

evaluated for data from Cyg X-l (dots) and Cyg X-3 (crosses) 
obtained in Goddard Flight 13.010 (see Figure 5) on 1973 
October 4. 


Figure 15 Principal frequencies (i.e. where a priori probability is 
£ 1C 4 for obtaining power density s observed value) below 
10 Hz for five consecutive 10.486 second intervals of data 
from Cyg X-l (Flight 13.010). Right margin presents a 
histogram of these principal frequencies (scale: 0-5). 
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Figure 


Figure 


Figure 


Figure 


16 Top: Principal frequencies (i.e., where a priori probability 

is £ 1CT 3 for obtaining power density observed value) 
above 10 Hz for consecutive 0.3277 second intervals of data 
from Cyg X-l (Flight 13.010). Right margin presents a his- 
togram of principal frequencies in the band 100 - 1500 Hz. 
Bottom: Counts per 320 microsecond bin (averaged over 

1024 bins) versus time (0.2377 second units) corresponding 
to data intervals used for above power spectral analysis. 

17 The normalized autocorrelation function [p( T ) defined by 
Brinkman et al. (1974)] for the low-energy bin (2. 1-5.1 keV) , 
calculated from the data of a typical UilURU sighting of 

Cyg 1 on 1971 June 9. The delay time x is in units of 
0.192s. (Brinkman et al., 1974). 

18 The normalized autocorrelation function p(x) for the lew- 
energy bin (2. 1-5.1 keV) and high-energy bin (5.1-12 keV) 
obtained from the summation of the autocorrelation functions 
for six UHURU sightings of Cyg X-l. (Brinkman et al., 1974). 

19 The normalized autocorrelation function p(r)» as defined 
by Brinkman et al. (1974), calculated from the Cyg X-l 
data of Goddard Flight 13.010 on 1973 October 4, taken 
separately for the energy bins 2. 1-5.1 keV and 5.1-12 keV. 

The delay time x *- s i- n units of 0.250s. 


32 


Figure 20 


Figure 21 


Figure 22 


Figure 23 
(top) 


Figure 23 
(bottom) 


The normalized autocorrelation function p ( t) calculated 
from the Cyg X-3 data of Goddard Flight 13.010, taken 
separately for the energy bins 2. 1-5.1 keV and 5.1-12 keV. 
The delay time t is in units of 0.205s. 

The normalized autocorrelation function p(x) calculated 
from computer generated shot-noise data based upon X = 20 
pulses/second of 0.5s duration each, noimialized to the 
average count rate of Cyg X-l during the exposure of Goddard 
Flight 13.010. The delay time \ is in units of 0.205s. 

The normalized autocorrelation function p(x) evaluated with 
0.0205s resolution for the Cyg X-l data from Goddard Flight 
13.010 and for three sets of computer generated shot-noise 
data based upon pulse widths of 0,5s, 0.1s, and 0.98s 
respectively. 

The normalized autocorrelation function pOr) calculated 
from the Cyg X-l data of Goddard Flight 13.010, taken 
separately for the energy bins 2. 1-5.1 keV, 5.1-12 keV and 
12-35 keV, with a resolution of 0.082s for the delay time. 

The normalized autocorrelation function p ( T ) calculated 
from the Cyg X-3 data of Goddard Flight 13.010 for the 
energy bin 2. 1-5.1 keV, with a resolution of 0„082s for 
the delay time. 
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Figure 24 A histogram of the observed population for the number of 

counts per 1.28 ms bin during the 409.6 ms duration interval 
#106 encompassing the peak Cyg X-l flux during Goddard 
Flight 13.010. The crosses give the expected Poisson 
distribution function corresponding to the average rate 
,(2.269 counts per 1.28 ms bin) and the associated statistical 
rms deviation. 

Figure 25 A histogram of the observed population for the number of 

counts per 1,28 bin during interval #106 of Goddard Flight 

13.010, but excluding the 20.48 ms duration frame #7. The 
crosses give the expected Poisson distribution function 
corresponding to the average rate (2.082 counts per 1.28 
bin) and the associated statistical rms deviation. 

Figure l6 Count rate profile (counts per temporal bin) for Cyg X-l 
during TM frames #6-8 of interval #106, Goddard Flight 

13.010. The telemetry word length of 320 jus is used as the 
finest resolution temporal bin, and the corresponding data 
are shown as the shaded histogram for frame #7. The over- 
all histogram shown for frames #6-8 is based on a 1.28 ms 


bin. 
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